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FOREWORD 


The  influence  of  finite  wall  impedance  effects  on  the  cyclotron  maser 
instability  of  a  hollow  electron  beam  is  investigated.  The  stability  analysis  is 
carried  out  within  the  framework  of  the  linearized  Vlasov-Maxwell  equations, 
assuming  that  the  beam  thickness  is  much  less  than  the  radius  of  the  beam.  The 
formal  dispersion  relation  for  azimuthally  symmetric  electromagnetic  perturbations 
is  obtained,  including  the  influence  of  an  arbitrary  value  of  wall  impedance. 

One  of  the  most  important  features  of  this  analysis  is  that,  for  a  purely 
resistive  wall,  the  growth  rate  of  instability  is  substantially  reduced  by 
allowing  even  a  very  small  amount  of  resistivity.  Moreover,  the  range  of  axial 
wave  numbers  corresponding  to  instability  increases  rapidly  as  the  wall  resis¬ 
tivity  is  increased.  Cyclotron  maser  stability  properties  in  a  dielectric 
loaded  waveguide  is  also  investigated.  It  is  shown  that  by  an  appropriate 
choice  of  the  dielectric  constant  e  and  thickness  of  dielectric  material,  the 
bandwidth  of  instability  can  be  increased  more  than  twicte  of  that  for  a  plain 
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I .  INTRODUCTION 


In  recent  years,  there  has  been  increasing  interest  in  the  electron 

1-9 

cyclotron  maser  instability  in  connection  with  intense  microwave 
generation.^  ^  For  the  most  part,  previous  theoretical  analyses  of 
this  instability  have  been  carried  out  without  including  the  influence  of 
finite  impedance  effects  of  a  waveguide  which  has  been  previously 
assumed  to  be  a  perfect  conductor.  Although  this  is  a  reasonable  assump¬ 
tion  in  the  present  experiments,  we  expect  significant  modification  to 
the  stability  behavior  when  a  small  amount  of  an  artificial  resistivity 
is  introduced  into  the  waveguide  wall  in  order  for  a  stable  operation  of 
the  microwave  amplification.  Moreover,  to  increase  the  bandwidth  of 
microwave  amplification,  a  dielectric  loaded  waveguide  is  more  desirable 
than  a  pure  conducting  waveguide.  However,  in  the  dielectric  loaded 
waveguide,  the  impedance  of  wall  is  purely  reactive.  In  this  paper,  we 
investigate  the  influence  of  finite  wall  impedance  effects  on  the  cyclo¬ 
tron  maser  instability  of  a  hollow  electron  beam  in  a  waveguide  with  an 
arbitrary  impedance  Z. 

Equilibrium  and  stability  properties  are  calculated  for  the  specific 
choice  of  electron  distribution  function  [Eq.  (3) ] . 

f°(H,P.,P  )  =  (to  N  /4TT2mc2)6(Y-f)6(p  -p)<5(P.-P  )  , 

eozce  z  oo 

where  H=ymc2  is  the  energy,  P.  is  the  canonical  angular  momentum,  P  is 

u  z 

the  axial  momentum,  to  =eB  /ymc  is  the  electron  cyclotron  frequency,  and 

c  o 

Ng  is  the  number  of  electrons  per  unit  axial  length.  The  stability  ana¬ 
lysis  in  this  paper  is  carried  out  within  the  framework  of  the  linearized 
Vlasov-Maxwell  equations,  assuming  that  the  beam  thickness  is  much  less 
than  the  equilibrium  radius  Rq  of  the  beam  and  that  v/^<<1,  where  v  is 
Budker's  parameter. 


7 


NSWC  TR  80-131 


The  formal  stability  analysis  for  azimuthally  symmetric  electromag- 
netric  perturbations  (3/39=0)  is  carried  out  in  Sec.  II,  including  the 
important  influence  of  finite  wall  impedance  Z,  which  in  general  is  an 
arbitrary  function  of  the  eigenf requency  uj  and  axial  wavenumber  k .  The 
dispersion  relation  of  cyclotron  maser  instability  in  Eq.  ( 20 i ,  when  com¬ 
bined  with  Eqs.  (17)  and  (18),  constitutes  one  of  the  main  results  of 
this  paper  and  can  be  used  to  investigate  stability  properties  for  a 
broad  range  of  system  parameters.  In  this  regard,  we  emphasize  that  Eq. 
(20)  is  derived  with  no  a  priori  assumption  that  the  electron  beam  is 
tenuous,  or  that  the  impedance  of  wall  is  independent  of  j  and  k 

In  Sec.  Ill,  sum  of  magnetic  wave  admittances  b_  +  b+  in  Eq.  (20) 
can  be  significantly  simplified  by  making  use  of  tenuous  beam  limit 
(v/y’+O) .  In  the  absence  of  a  oeam,  the  vacuum  transverse  electric  (TE) 
mode  dispersion  relation  is  given  by  Eq.  (26) 

,  2  .  2  ,  2.  2  2  ,  ,  N 
(ui  /c  -k  )R  =x  (uj,k)  , 
w  on 

where  c  is  the  speed  of  light  in  vacuo ,  R^  is  the  radius  of  impedance 

wall,  x  is  the  nth  root  of  Z=J, (x  )/x  J  (x  ) ,  and  J  (x)  is  the 
on  1  on  on  o  on 

Bessel  function  of  the  first  kind  of  order  i.  The  root  x  is  calculated 

on 

numerically  for  a  given  complex  value  of  wall  impedance  Z.  The  detailed 
dependence  of  the  root  x  in  Eq.  (26)  on  the  eigenfrequency  oi  and  axial 
wavenumber  k  is  particularly  important  in  connection  with  the  gain  and 
bandwidth  of  a  microwave  amplification  by  the  cyclotron  maser  instability. 

In  Sec.  IV,  a  detailed  analytic  and  numerical  investigation  of  the 
cyclotron  maser  instability  is  carried  out  for  a  wall  impedance  indepen¬ 
dent  of  the  eigenfrequency  w  and  axial  wavenumber  k  (i.e.  ,  3Z/3w=0, 
3Z/3k=0)  .  Defining  the  wall  impedance  Z  by  Z=  |  Z  |  exp  (-i<f> )  ,  stability  pro¬ 
perties  of  cyclotron  maser  instability  are  investigated  for  a  broad  range 
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of  the  magnitude  jz,  and  the  phase  angle  $ .  One  of  the  most  important 
features  of  this  analysis  is  that,  for  a  purely  resistive  wall  charac¬ 
terized  by  $=90°,  the  growth  rate  is  substantially  reduced  by  introducing 
even  a  very  small  amount  of  resistivity  into  the  wall.  This  feature 
represents  a  general  tendency  for  all  radial  mode  number  n.  Moreover, 
the  range  of  the  normalized  axial  wavenumber  kc/j^  corresponding  to 
instability  increases  rapidly  as  the  resistivity  of  wall  is  increased. 

Cyclotron  maser  stability  properties  in  a  dielectric  loaded  wave¬ 
guide  are  investigated  in  Sec.  V,  assuming  that  the  impedance  of  dielec¬ 
tric  material  is  purely  reactive.  It  is  shown  that  stability  properties 
of  cyclotron  maser  instability  exhibit  a  sensitive  dependence  or.  the 
dielectric  constant  e  and  thickness  of  dielectric  material.  By  an  appro¬ 
priate  choice  of  e  and  thickness  of  dielectric  material,  the  bandwidth  of 
instability  can  be  increased  more  than  twice  of  that  for  a  plain  con¬ 
ducting  waveguide.  Moreover,  dielectric  material  reduces  the  radial 
wavelength  of  the  waveguide  mode,  thereby  significantly  modifying  the 
nature  of  perturbations.  The  influence  of  self-dielectric  effects  of  an 
electron  beam  itself  on  the  stability  behavior  is  investigated.  However, 
for  a  low  density  beam  (i.e.,  v/f<<0.005)  with  small  transverse  velocity, 
self-dielectric  effects  are  insignificant. 
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LINEARIZED  VLASOV-MAXWELL  EQUATIONS 


The  equilibrium  configuration  consists  of  a  hollow  relativistic 
electron  beam  propagating  parallel  to  a  strong,  externally  applied  mag¬ 
netic  field  B„e  .  The  radius  of  the  electron  beam  is  denoted  by  R  ,  and 
0-Z  o 

a  finite  impedance  wall  is  located  at  radius  r=R  .  Cylindrical  polar 

w 

coordinates  (r,6,z)  are  introduced.  In  the  present  analysis,  we  assume 
that 

*/'?"<  1  ,  (1) 

2  2 

where  v=N  e  /me  is  Budker 1 s  parameter, 
e 

R  o 

N  =2tt f  edrr  n  (r)  ,  (2) 

e  o  e 

is  the  number  of  electrons  per  unit  axial  length,  n°(r)  is  the  equilibrium 

electron  density,  c  is  the  speed  of  light  in  vacuo ,  -e  and  m  are  the 

2 

electron  charge  and  rest  mass,  respectively,  and  ?mc  is  the  characteris¬ 
tic  electron  energy  in  the  laboratory  frame.  Consistent  with  the  low- 
intensity  assumption  in  Eq.  (1),  we  also  neglect  the  influence  of  equili¬ 
brium  self-field. 

In  the  present  article,  we  investigate  the  stability  properties 
associated  with  the  beam  distribution  function 

w  N 

f°(H,Pg,Pz)=  — S(Y-?)6(pz-p)6(P9-Po)  ,  (3) 

4ti  ■? me 

2  2422  1/2 

where  H=ymc  = (m  c  +c  p  )  is  the  total  energy,  PZ=PZ  ls  the  axial 

canonical  momentum,  P  =r[P  -(e/2c)rB  ]  is  the  canonical  angular  momentum, 

0  0  o 

u)  =  eB  /me  is  the  nonrelativistic  electron  cyclotron  frequency, 
c  o 

P  =  - (e/2c) (R2-r^)B  ,  (4) 

o  o  L  o 

is  the  canonical  angular  momentum  of  an  electron  with  Larmor  radius, 

2  2  "2  '  "  2  1/2 

r  =  [  (y  -l)c  /u>  -  (p/mu)  )  \  '  ,  (5) 

L  c  C 
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and  p  is  a  constant.  For  a  notational  convenience  in  the  subsequent 

2  2  1/2 

analysis,  we  also  introduce  the  definitions  8  =p/ymc  and  6, =(1-8  -1/y  ) 

z  z 

The  equilibrium  properties  of  an  electron  beam  described  by  Eq.  (3)  can 
be  found  in  the  previous  literatures. 3 ' ^ 

In  the  subsequent  stability  analysis,  we  make  use  of  the  linearized 
Vlasov -Maxwell  equations  of  azimuthally  symmetric  perturbations  0/30=0) 
about  a  tenuous,  hollow  beam  equilibrium  described  by  Eq.  (3).  We  adopt 
a  normal-mode  approach  in  which  all  perturbations  are  assumed  to  vary 
with  time  and  z  according  to 

6iMx,t)=<Mr)exp[i  (kz-wt)  ]  , 

where  Imw^O  Here,  w  is  the  complex  eigenf requency  and  k  is  the  axial 
wavenumber.  The  Maxwell  equations  for  the  perturbed  electric  and  magne¬ 
tic  field  amplitudes  can  be  expressed  as 

V  x  E(x)=i(w/c)B(x)  , 

V  x  (l/u)B  (x)=  (4ir/c)  J  (x) -i  (u/c)  e  E(x)  ,  (6) 

where  e  and  u  are  the  dielectric  constant  and  permeability,  respectively, 
of  the  background  material,  E(x)  and  B(x)  are  the  perturbed  electric  and 
magnetic  fields,  respectively,  and 

J(x)=  -e /d3p  v  f  (x,p)  (7) 

is  the  perturbed  current  density.  Note  that  E=y=l  in  vacuo.  In  Eq.  (7)  , 

f  (x,p)=e f°  dx  exp(-iwT) 

0  —  CO  — 

is  the  perturbed  distribution  function,  T=t'-t,  and  the  particle  trajec¬ 
tories  x'(t')  and  p'(t')  satisfy  dx/dt'=v'  and  dp'/dt'=  -ev'  x  b  e  /c, 

-0_Z 

with  "initial"  conditions  x'(t'=t)=x  and  v'(t'=t)=v. 


E  (x')+  ~ 


x  B(x') 


(8) 


3P' 
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Making  use  of  Eq.  (6),  it  is  straightforward  to  show  that 


B  (r)  =  -  (kc/w) Efj  {r )  , 


B  (r)=  -i  (c/ur )  (  d  (rE., ) /lr  )  , 

z  0 


—  fi.  3-  (rt  )1  +  p2e 

<r  L  ur  dr  6  J  t 


4TTjg)p  ‘ 

c  0 


2  2.2,2 

p  =«  in  /c  -k  ,  (11) 

!■:  is  the  azimuthal  component  of  perturbed  electric  field,  and  B  and  B 

r 

are  the  radial  and  axial  components,  respectively ,  of  the  perturbed  mag¬ 
netic  field.  The  field  equations  (9)  and  (10)  represent  the  transverse 

electric  (TE)  waveguide  modes.  For  the  cyclotron  maser  instability,  it 
1  9 

is  well  established  '  that  the  beam  cyclotron  resonance  mode  couples 
much  more  strongly  with  the  TE  waveguide  mode  than  the  transverse  magne¬ 
tic  (TM)  waveguide  mode.  In  this  regard,  the  present  analysis  is  restric¬ 
ted  to  the  TE  waveguide  mode . 

The  perturbed  azimuthal  electric  field  E^tr)  is  continuous  across 
the  beam  inner  and  outer  boundaries  r=R^  and  r=R2 ’  Integrating  Eq.  (10) 
from  r=R^-6  to  ^R^+fi  and  taking  the  limit  6-*-0+,  we  obtain  the  approxi¬ 


mate  result 


+ 

R 

Bz(R2)_Bz(Ri)=  ~  H  -  drJ0(r)  ' 

R1 


where  iHftr-)  denotes  lim,  „  i(i(R_+6),  and  use  has  been  made  of  the  hollow 
2  o-K)  2— 

+ 

beam  approximation  (Thickness  of  the  beam  is  much  less  than  the  beam 
radius  Rq) .  For  convenience  in  the  subsequent  analysis,  we  introduce 
the  normalized  magnetic  wave  admittance^-  b+  defined  at  the  inner  and 
outer  surfaces  of  the  electron  beam  by 
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b  =  -iu/p.cR  ) B  (R  ) /E  (R  )  , 

-  i  O  Z  I  O 


b  =i(w/p.cR  )  B  ( R  )  /E  (R  )  , 

+  1  o  z  2  6  o 


(13) 


where 


2  2.2,2 

PL=^  /c  -k  . 


(14) 


Moreover,  we  also  define  the  wave  impedance  Z(u>,k)  of  the  wall  as 


ZU,k)=  -i  (c/wR  )  E,  (R  )/B  (R  )  , 

w  -J  w  z  w 


(15) 


where  R^  is  the  radius  of  an  impedance  wall. 

Since  the  perturbed  current  density  vanishes  in  the  vacuum  region 
outside  the  beam,  it  is  obvious  that  the  solution  to  Eq.  (10)  in  the 
vacuum  region  is  given  by 


Vr): 


A  J1 (p^r ) ,  r<R1  , 


B  J^p^J+C  N1(p1r),  R2<r<Rw, 


(16) 


where  J.(p.r)  and  N. (p,r)  are  Bessel  functions  of  the  first  and  second 

l  *i  V.  1 

kind,  respectively,  of  order  l.  Making  use  of  Eqs.  (15)  and  (16) ,  and 
the  thin  beam  approximation  (R^-^i^^iO  '  we  exPress  the  sum  k_+k+  °f  the 
wave  admittance  in  Eq.  (13)  as 


b  +b  = 


[2g(U  A)  (R  /R  C)  ‘ 
w  o 


+  J  (CR  /R, ,)  [j,  (CR  /R  )+g(C)N  (CR  /R  )  ] 
lowiow  low 


(17) 


where  C=P.R  ,  and 
1  w 


g(C)  = 


CJ  (C) 

o 


n^O-Czn^C) 


J1(5) 

TTJt) 


(18) 


Substituting  Eq.  (13)  into  Eq.  (12),  we  obtain 

2 


4ttwR 


<b-+V=  -1  ~ 2T2 


w 


c  C  R  E  (R  ) 
o  0  o 


dr  J0(r) 


(19) 
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After  some  tedious  but  straightforward  algebraic  manipulation  with  Eqr . 
(7) ,  (8)  and  (9) ,  we  obtain  the  dispersion  relation  of  the  cyclotron 
maser  instability 

2 


2,22 
a)  -k  c 

(w-kB  c-u  ) 
z  c 


2  ' 


(20) 


where  u  =w  /y  is  the  electron  cyclotron  frequency  on  the  laboratory  frame, 
c  c 


and  use  has  been  made  of  the  assumption  6A<<1.  For  a  detailed  deriva¬ 


tion  of  Eq.  (20) ,  we  urge  the  reader  to  review  the  previous  literatures. 
In  the  remainder  of  this  article,  we  investigate  the  stability  properties 
of  Eq.  (20)  for  various  wall  impedance  Z  on  Eq.  (18) . 


1,3,5 


*~r 


TT 
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r 1 1 .  LIMIT  OF  A  TENUOUS  BEAM 


To  simplify  the  expressions  for  the  sum  b  +b+  of  the  wave  admittance 
in  Eq.  (17),  it  is  useful  to  consider  Die  limit  where  the  beam  density 
is  very  tenuous,  i.e., 

v/- y'-O  (21) 


In  this  absence  of  a  beam,  the  vacuum  TE:  mode  dispersion  relation  is 
UlVi'ti  by 

Z=J.(U/»,J  ((,)  ,  (22) 

I  o 

where  Z(to,k)  is  the  wave  impedance  of  the  wall  [Eq.  (15))  and 
222  1/2 

£=(to  /c  -k  )  R^.  Substituting  Eq.  (22)  into  Eq.  (18),  we  can  simplify 

Eq.  (17)  by 


b  +b  =  - 
-  + 


r  r  j  (o  i 

2 

r 

w  o 

i 

R  Jn (£R  /R  ) 

Z  f,  J  ( r ) 

o  1  o  w 

L  °  J 

(23) 


For  a  convenience  in  future  analysis,  we  define  the  nth  root  of  Eq.  (22) 
by 


£=x  (<o,k)=x  (w,k)-i  x. (w,k)  (24) 

on  r  i 

for  a  specified  wall  impedance 

Z=|  z|exp(-i<)i)=Zr+iZi  .  (25) 

For  a  very  tenuous  beam  (Eq.  (21) J,  it  is  evident  from  Eqs.  (22)  and 
(24)  that 


(to2/c2-k2)  R2=x2  (w,k)  (26) 

w  on 

is  a  good  approximation  to  the  dispersion  relation  in  Eq.  (20) .  The 
right-hand  side  of  Eq.  (20)  describes  the  beam-produced  modifications 
to  the  vacuum  dispersion  relation  (26).  We  remind  the  reader  that  the 
nth  root  x  is  a  function  of  the  eigenfrequency  to  and  axial  wavenumber 
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k  in  general.  For  an  impedance  with  a  small  magnitude  satisfying  |Z|<<1, 
we  obtain 

x  =a  [ 1+ | Z | exp (-i$) ]  (27) 

on  on 

from  Eqs.  (22)  and  (25).  Therefore,  the  vacuum  waveguide  mode  for  zero 
wall  impedance  can  be  expressed  as 

2  2  2  2  2 

(u>  /c  -k  )R  =cT  (28) 

w  on 

where  a  is  the  nth  root  of  J, (a  )=0,  which  is  independent  of  uj  and  k. 
on  1  on 

Figure  1  is  the  contours  of  constant  phase  angle  $  and  magnitude 

I Z I  of  the  function  Z=J, (x  ) /x  J  (x  )  in  the  plane  of  the  root  x  = 

1  on  on  o  on  on 

x^-ix^  for  (a)  n=l  and  (b)  n=2.  We  note  from  Fig.  1  that  the  nth  root 

x^^  of  Eq.  (22)  approaches  to  a  and  the  magnitude  of  wall  impedance  is 

reduced  to  zero.  On  the  other  hand,  the  nth  root  x  approaches  to  6 

on  on 

by  increasing  the  magnitude  |z|  to  infinity.  Here  6  is  the  nth  root 

of  J  (8  )  =0.  We  also  emphasize  that  for  specified  values  of  Izl  and  d> , 

o  on 

the  root  x  =x  -ix.  is  determined  from  Fig.  1.  For  example,  for  |z|=0.3, 
on  r  i  11 

<t>=90°  and  n=l,  we  find  from  Fig.  1  that  xr=3.12  and  x^=0.82.  Obviously, 
the  root  XQn  of  Eq.  (22)  is  a  very  complicated  function  of  the  wall  impe¬ 
dance  Z. 

The  dependence  of  the  nth  root  xq^  on  the  eigenfrequency  io  and  axial 
wavenumber  k  is  particularly  important  in  connection  with  the  gain  and 


bandwidth  of  a  microwave  amplification  by  the  cyclotron  maser  instability. 

2  2  2  2  2  1/2 

Shown  in  Fig.  2  are  schematic  drawings  of  the  curve  io= (k  c  +a  c  /R  ) 

’  on  w 

2  2  2  2  2  1/2 

corresponding  to  zero  wall  impedance  and  the  curve  w=(k  c  +xQnc  /R  ) 
corresponding  to  an  arbitrary  wall  impedance.  The  straightline  u)=k6zc+(oc 
is  the  cyclotron  resonance  mode.  It  is  important  to  note  from  Fig.  2 


that  the  interaction  region  of  cyclotron  resonance  mode  with  vacuum  wave¬ 


guide  dispersion  relation  is  increased  substantially  by  an  appropriate 
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choice  of  the  wall  impedance  (x  ) .  In  this  regard,  we  conclude  that 
the  growth  rate  or  the  bandwidth  of  cyclotron  maser  instability  can  be 
considerably  improved  when  the  functional  dependence  of  the  wall  impe¬ 
dance  on  u  and  k  is  matched  with  the  cyclotron  resonance  mode  in  a  broad 

range  of  k  space.  However,  a  detail  investigation  of  wall  impedance  is 

14 

required.  Previous  studies  indicate  that  a  proper  waveguide  structure 
provides  a  necessary  wall  impedance.  The  influence  of  nonuniform  wall 
impedance  on  cyclotron  maser  instability  is  currently  under  investigation 
by  authors  for  a  broad  range  of  physical  parameters  and  waveguide  struc¬ 
tures.  As  a  simple  example,  we  present,  in  Sec.  V,  the  cyclotron  maser 
instability  in  a  dielectric  loaded  waveguide. 


19/20 
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I V .  CYCLOTRON  MASER  INSTABILITY  IN  A  WAVEGUIDE  WITH 
CONSTANT  WALL  IMPEDANCE 


The  cyclotron  maser  stability  properties  of  a  hollow  electron  beam 

in  a  vacuum  waveguide  with  general  impedance  are  investigated  in  this 

section,  assuming  that  impedances  of  the  wall  are  independent  of  the 

eigenfrequency  u>  and  the  axial  wavenumber  k,  i.e., 

3Z/3w=0,  3z/3k=0  (29) 

2  2  2  2  2 

Taylor  expanding  Eq.  (23)  about  (u,  /c  -k  ) Rw=xQn  with  Eq.  (29),  it  is 
straightforward  to  show  that  b  +b+  can  be  approximated  in  leading  order 
by 


b  +b  = 
—  + 


2x 


on 


J  (x  )  -J  (x  )  J  (x  ) 
1  on  o  on  2  on 

(x  R  /R  ) 

1  on  o  w 


(30) 


Substituting  Eq.  (30)  into  Eq.  (20)  and  making  use  of  Eq.  (24)  ,  we  obtain 


the  approximate  dispersion  relation 


2  „  x 

u  2_  on 

2  "  ~  2 

c  R 

w 


vB. 


J.  ( X  R  /R  ) 
1  on  o  w 


2.22 
oo  -k  c 


?R2  J2  (x  )-J  (x  )  J.  (x  )  (w-kB  c-oo  )2 
w  1  on  o  on  2  on  z  c 


(31) 


for  a  tenuous  thin  beam.  Stability  analysis  of  the  dispersion  relation 
in  Eq.  (31)  is  carried  out  in  the  remainder  of  this  section  for  a  broad 
range  of  physical  parameters. 

In  order  to  investigate  the  influence  of  wall  resistivity  on  stability 

behavior,  the  analysis  of  the  dispersion  relation  in  Eq.  (31)  is  carried 

out  for  a  wall  impedance  with  a  small  magnitude  satisfying  | Z | < < 1 .  In 

this  limit,  we  choose  the  wall  radius  R^  according  to 

R  =a  c/w  y  (32) 

w  on  c  z 

2  2-1 

for  the  present  purposes.  Here  Y  is  defined  by  Y  =(1-8  )  .  Equation 


(32)  ensures  that  the  group  velocity  of  the  vacuum  waveguide  mode  for 
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| Z J  =0  is  equal  to  the  beam  velocity  for  a  certain  k  value. ^  Making  use 

of  Eq.  (32)  and  defining  the  normalized  Doppler-shifted  eigenfrequency 

X  =  (oj-k8  c-ui  )/u  ,  (33) 

zee 

it  is  straightforward  to  show  that  the  dispersion  relation  in  Eq.  (31) 
can  be  expressed  as 


The  normalized  growth  rate  /.=lmx  and  Doppler- shifted  real  frequency 

Xr=Rex  have  bven  calrul ited  numerically  from  Eq.  (34)  for  v=0.002,  y=1.118, 

Bjl=0.4  and  several  different  values  of  wall  impedance  Z.  For  specified 

values  of  Z,  the  root  x  is  determined  from  Fig.  1.  For  example,  we 

on 

find  from  Fig.  1  that  x  =(x  ,x.)  =  (3.82,  0.167)  for  Z=  (  |  Z  |  ,  <J>)  =  (0.045, 

on  r  l 

90°)  . 

Although  the  present  experiments  of  the  cyclotron  maser  instability 

are  carried  out  in  a  lossless  cylindrical  conductor,  it  is  necessary  to 

investigate  the  stability  properties  for  a  resistive  waveguide,  since  an 

artificial  resistivity  of  the  wall  is  sometimes  required  for  a  stable 

operation  of  the  microwave  amplification.  In  this  regard,  we  investigate 

the  dispersion  relation  in  Eq.  (34)  for  a  pure  resistive  wall  characterized 

by  4>=90°.  Shown  in  Fig.  3  are  plots  of  (a)  the  normalized  growth  rate 

and  (b)  the  normalized  Doppler-shifted  real  frequency  Xr  versus  kc/a)c 

obtained  from  Eq.  (34)  for  n=l,  R  /R  =0.5,  <}>=90o,  and  several  values  of 

o  w 

I Z | .  In  Fig.  3(b),  xr  is  plotted  only  for  the  ranges  of  kc/w,  corres¬ 
ponding  to  instability  (x^>0)  .  Several  features  are  noteworthy  from  Fig. 

3.  First,  the  maximum  growth  rate  decreases  as  the  magnitude  |z|  of  wall 
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impedance  is  increased.  The  growth  rate  is  substantially  reduced  by 
introducing  even  a  very  small  amount  of  resistivity  into  the  wall.  As 
evident  from  Eq.  (34)  and  Fig.  1,  this  feature  represents  a  general  ten¬ 
dency  for  all  radial  mode  numbers  n.  Second,  the  range  of  kc/i  corres- 

c 

ponding  to  instability  increases  rapidly  as  the  resistivity  of  the  wall 
is  increased.  This  mechanism  somehow  broadens  the  bandwidth  of  microwave 
amplification.  Third,  for  lossless  conductor,  the  Doppler-shifted  real 
frequency  increases  abruptly  when  the  wavenumber  k  approaches  to  the 
k-space  boundary  corresponding  to  marginal  stability  (.^=0).  However, 
resistivity  of  the  wall  levels  the  Doppler-shifted  real  frequency  [Fig. 
3(b) ] . 

The  dependence  of  stability  properties  on  phase  angle  <j>  of  the  impe¬ 
dance  is  illustrated  in  Fig.  4  where  the  normalized  growth  rate  is 

plotted  versus  kc/w  for  n=l,  R  /R  =0.5,  I Z 1=0. 05  and  several  different 
c  o  w 

phase  angles  <i>.  Since  the  plot  of  the  Doppler-shifted  real  frequency 
has  a  form  similar  to  that  of  Fig.  3(b),  it  is  not  shown  in  Fig.  4. 
Evidently  from  Fig.  4,  the  maximum  growth  rate  occurs  at  the  phase  angle 
4>=180°.  For  this  particular  magnitude  of  impedance  ( | Z [ =0 . 05 ) ,  the  growth 
rate  of  instability  is  reduced  sharply  as  the  phase  angle  of  impedance 
is  reduced  to  zero. 

To  investigate  the  cyclotron  maser  instability  for  a  purely  reactive 

impedance15  characterized  by  4>=0°  or  $=180°,  we  analyze  the  dispersion 

relation  in  Eq.  (31)  for  different  values  of  x  with  x.=0.  In  order 

on  i 

to  maximize  the  growth  rate  and  efficiency  of  microwave  amplification 
we  choose1 


V*o»c/“cV  VV-lF-on  • 


<351 


where  a  is  the  first  root  of  and  (x)=dJ^/dx.  Making  use 

of  Eqs.  (33)  and  (35),  we  can  show  that  the  dispersion  relation  in 
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Eq.  (31)  is  expressed  as 


Figure  5  presents  plots  of  the  normalized  maximum  growth  rate  (solid  line) 

and  R  /R  (broken  line)  versus  x  =x  obtained  from  Mas.  ( 36)  and  (  IS) , 
o  w  on  r 

respectively,  for  v=0.002,  y=1.118,  and  6^=0. 4.  Numerical  calculation 

indicates  that  the  maximum  growth  rate  occurs  at  the  wavenumber  k -0 . 15e  /c 

c 

for  the  range  2^x  ^_4  as  shown  in  Fig.  5.  We  note  from  Fig.  5  that  the 

growth  rate  of  instability  is  substantially  increased  by  reducing  the 

value  x  to  two.  Moreover,  it  may  be  possible  that  further  increase  of 
on 

the  growth  rate  can  be  attainable  by  reducing  the  value  xQn  less  than  two 
and  by  arranging  suitably  the  beam  radius  so  that  electrons  never  hit  the 
wall  (see  Eqs.  (31)  and  (36)]. 
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V .  CYCLOTRON  MASER  INSTABILITY  IN  A  HI ELECTRIC  LOADED  WAVE'JL’IOE 

In  this  section,  we  investigate  the  cyclotron  maser  stability  pro¬ 
perties  of  a  hollow  electron  beam  in  a  dielectric  loaded  waveguide.  As 
illustrated  in  Fig.  6,  a  hollow  electron  beam  with  radius  is  propaga¬ 
ting  through  a  cylindrical  dielectric  material  with  its  inner  radius  R  . 


A  grounded  cylindrical  conducting  wall  is  located  at  radius  R 


In 


general,  the  permeability  ;i  of  a  dielectric  material  differs  from  unity 

5 

by  only  a  few  parts  in  10  (u^l  for  paramagnetic  substances,  1  tor  dia¬ 
magnetic  substances) .  Therefore,  we  make  the  approximation  u=l  in  the 
calculation  in  this  section.  In  this  regard,  the  perturbed  azimuthal 

electric  field  (r)  and  axial  magnetic  field  B  (r)  are  continuous  across 
9  Z 

the  a  _4.ectric  boundary  r=R  .  From  Eq.  (10)  ,  we  obtain 

w 


3  f  i  3  / \ 
3r  [r  3r  (rLe) 


+p2e0=o  , 


inside  the  dielectric  material.  Here  the  parameter  p2  is  defined  by 

2  2.22 
P2=uj  e/c  -k 

The  solution  to  Eq.  (37)  can  be  expressed  as 

EQ  (r)=A(J1  (p2r)  -Jj  (C)N1  (p^)/^  (£)  )  , 


(37) 


(38) 


(39) 


2  22  1/2 

where  A  is  an  arbitrary  constant  and  C=RcP2=Rc (u>  e/c  -k  )  .  Obviously 

Eq.  (39)  satisfies  the  boundary  condition  Eg(r=Rc>=0.  Substituting  Eq. 
(39)  into  Eq.  (9)  and  making  use  of  the  boundary  conditions  at  r=Rw»  we 
obtain  the  wave  impedance  Z(ui,k) 


Z  ( oo ,  k. )  - 


J1(n)N1(c)-J1(ON1(n) 

n (J  (n)N  (;)-j.  (ON  (n)  1  ' 

o  1  1  o 


(40) 


from  the  definition  in  Eq.  (15).  In  Eq.  (40),  the  parameter  n  is  defined 
2  22  1/2 

by  n=RwP2=Rw((»)  e/c  -k  )  .  Note  from  Eq.  (40)  that  the  wall  impedance 

of  dielectric  loaded  waveguide  is  a  function  of  eigenfrequency  to  and 
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axial  wavenumber  k.  Since  the  dielectric  constant  e  is  always  greater 

than  unity,  we  note  from  Eqs.  (14)  and  (38)  that  p^-p^.  We  therefore 

conclude  from  the  definitions  t=p,R  ,  n=p„R  and  ;=p  R  that 

1  w  2  w  2c 

R  /R  =n/C  (41) 


C<n<_;  .  (42) 

Equation  (20),  combined  with  Eqs.  (17),  (18)  and  (40),  constitutes  a  com¬ 

plete  dispersion  relation  of  cyclotron  maser  instability  in  a  dielectric 
waveguide . 

Previous  study1  indicates  that  the  maximum  growth  rate  occurs  at 

the  sum  of  wave  admittance  b_+b+=0  characterized  by  the  vacuum  TE  mode 

dispersion  relation  [Eq,  (22)]  in  the  region  0<r<Rw-  Making  use  of  the 

definitions  5=p  R  ,  n=p^R  and  £=p  R  ,  we  can  express  the  vacuum  TE  mode 
1  w  2  w  2  c 

dispersion  relation  on  Eq.  (22)  as 

2222  222  22 

u/c-k-x*  /R  =r/R  =(€C/n)  /R  ,  (43) 

on  w  w  c 

where  the  parameter  £  is  related  to  n  by 

J  (n)N  (?)-J  (ON  (n) 

Z=J  (O/CJ  (£)  = - - -  ,  (44) 

°  n(J  (n)N  (c)-j.  (ON  (n)l 

o  1  1  o 

for  a  specified  value  of  parameter  £.  It  is  instructive  to  examine 

Eq.  (43)  in  the  limit  e-*-l.  Making  use  of  £=h>  and  Eq.  (44),  we  have 

(oj2/c2-k2)  R2=a2  from  Eq.  (43)  for  e=l. 
c  on 

Figure  7  is  plots  of  the  impedance  Z  versus  parameter  £  (broken  lines) 
and  Z  versus  parameter  n  (solid  lines)  obtained  from  Eq.  (44)  for  £=4,  5 
and  6.  The  scales  in  horizontal  line  in  Fig.  7  represent  both  the  para¬ 
meters  £  and  n.  The  solid  curves  are  plotted  only  for  the  ranges  n  satis¬ 
fying  Eq.  (42).  Note  that  for  specified  values  of  £  and  C,  the  parameter 
h  is  determined  from  Fig.  7.  For  example,  for  £=2,  we  find  from  Fig.  7 
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that  7=2.175  for  ',  =  4,  n=3.3  for  ;=5  and  ”  =  4.35  for  '=6  i-orri'^or.dir...:  to 
the  wave  impedance  Z=1.25  at  surface  of  a  dielectric  material.  1::  this 
particular  example  (:,=2)  ,  values  of  parameters  -l:,/”  and  R  K  are  riven 

W  C.' 

by  (£C/n,  R  /R  )  =  {3.68,  0.54)  for  ".=4,  (3.03,  ).66)  for  1  =  5  and  (2.70, 

w  c 

0.73)  for  1=6.  From  the  definitions  of  1,  ”  and  1,  we  obtain  .  (  -1)= 

2  2  2  2  2 

c  (!-•"  /n  )1  /R  which  unmistakably  indicates  that  the  eigenf requency  of 
vacuum  waveguide  mode  depends  on  the  dielectric  constant  t.  However, 

the  maximum  growth  rate  of  cyclotron  maser  instability  of  electron  beams 

2  2  ■>  2 

characterized  by  ^“-•'•'1  occurs  at  the  eigenf  requency  -  satisfying  , 

2  2  2  2 

thereby  approximating  n  /'  =p^/p  by 


e=n2/^2  .  (45) 

In  this  regard,  we  can  approximately  determine  the  dielectric  constant 

e  from  Fig.  7  and  Eq.  (45)  for  the  cyclotron  maser  instability.  As  an 

example,  for  C=2,  we  obtain  e=1.19  for  £=4,  e=2.62  for  C=5,  and  e=4.75 

for  C=6.  We  also  emphasize  that  the  impedance  Z  on  Eq.  (40)  is  purely 

reactive  [i.e.,  Z^O  where  Z^  is  defined  in  Eq.  (25)]  for  the  eigenfre- 

2  2  2 

quency  ui  satisfying  u  >k  c  . 

In  order  to  complete  the  stability  analysis  of  cyclotron  maser 
instability  on  a  dielectric  loaded  waveguide,  it  is  required  to  investi¬ 
gate  numerically  the  dispersion  relation  in  Eq.  (20) ,  where  no  a  priori 
assumption  is  made  that  the  beam  is  very  tenuous.  However,  use  is  made 
of  the  fact  that  the  magnitude  of  Doppler-shifted  eigen frequency  |oo-kBzc- 

u  is  much  less  than  the  electron  cyclotron  frequency  (i.e.,  w-kB  c-u>  <<io  ). 
c  zee 

Evaluating  the  parameters  C»  h,  and  the  wave  admittance  b  +b+  at  ui=ujo= 

kB  c+u>  ,  the  dispersion  relation  in  Eq.  (20)  can  be  approximated  by 
z  c 


r  d 

1  3  ^  .  2  v 

6i2  /* 

1  \2 

W  l 

\l  \2  A 

u)  —  (b  +b  ) 
1  c  dai  -  + 

|  Y  +  (b  +b  )  X  +  fTT 
w  *  -  +  to  *  2-9 

Jo  o 

-J 

1  1+  8  z  kc/u>  1  -(kc/ui^) 

=0  (46) 
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where  '^=\  (~=~  ),  and  the  normalized  Doppler-shifted  eigenf requency  ;<  is 

defined  in  Eq.  (33).  The  normalized  growth  rate  x^=lmx  has  been  calculated 

numerically  from  Eq.  (46)  for  a  broad  range  of  system  parameters  kc/iii  ,  R  u  /c, 

c  w  c 

R  -^/c  and  dielectric  constant  c.  In  this  numerical  calculation,  it  has 

been  assumed  that  the  beam  parameters  are  v=0.002,  f=1.118  and  8^=0. 4. 

To  maximize  the  growth  rate  of  instability,  we  also  assumed^  R  w  /c=a, , , 

o  c  11 

where  a  is  the  first  root  of  J  1  (a  )=0. 

Fig.  8  illustrates  influence  of  the  self-dielectric  effects  of  an 

electron  beam  on  the  cyclotron  maser  instability.  Shown  in  Fig.  8  are 

plots  of  normalized  growth  rate  versus  kc/u ^  obtained  from  Eq.  (34) 

(broken  line)  for  x  =u„,  and  x  R  /R  =a, , ,  and  from  Eq.  (46)  (solid 
on  01  on  o  w  11 

line)  for  t  =  l,  R  ui  /c=R  u)  /c=a„  /y  =3.75.  Obviously  from  Fig.  8,  the 
wcccOl  z  1  ” 

growth  rate  calculated  from  the  self-consistent  dispersion  relation  in 

Eq.  (46) ,  including  the  self-dielectric  effects  of  beam  itself,  is 

slightly  greater  than  the  growth  rate  from  the  approximate  dispersion 

relation  in  Eq.  (34).  However,  for  a  low  density  beam  with  small  trans- 
2 

verse  velocity  (8^  <<1)  as  shown  in  Fig.  8,  self-dielectric  effects  of 
beam  itself  is  insignificant. 

Stability  boundaries  in  the  parameter  space  (R  w  /c,  c)  are  illus- 

c  c 

trated  in  Fig.  9  for  several  values  of  R^w^/c  ranging  from  2.2  to  2.8. 

In  Fig.  9,  the  solid  curves  correspond  to  the  stability  boundaries 
obtained  from  Eq.  (46)  for  ‘9=1.118,  8  =0.4.  For  a  given  value  of  R  w  /c, 

A  W  C 

the  region  of  (R  u  /c,  e)  parameter  space  below  the  curve  corresponds  to 
c  c  - 

absolute  stability  for  any  arbitrary  values  of  normalized  axial  wavenumber 

kc/w  ,  whereas  the  reqion  of  paremeter  space  above  the  curve  corresponds 
c  - 

to  instability  for  some  values  of  wavenumber  kc/u^.  Several  points  are 
noteworthy  from  Fig.  9.  First,  for  general  values  of  the  parameter  R^w^/c, 
the  stability  boundaries  converge  to  the  common  value  Rcwc/c~3.6  for  e=l 
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corresponding  to  vacuum  dielectric  constant.  In  fact,  the  stability 

properties  of  Eq.  (46)  are  the  same  for  any  arbitrary  values  of  R^uj^/c 

for  e=l.  Second,  for  a  given  radius  of  dielectric  material,  the  value  of 

dielectric  constant  e  required  for  instability  increases  rapidly  to  infinity 

as  the  radius  of  conducting  wall  (R^a  /c)  approaches  to  the  inner  radius 

of  dielectric  material  (R  *^/c).  Third,  obviously  for  a  given  radius  of 

conducting  wall,  instability  occurs  even  at  small  values  of  dielectric 

constant  as  the  value  of  parameter  R  ^  /c  is  decreased. 

w  c 

The  dependence  of  stability  properties  on  dielectric  constant  e  is 

further  illustrated  in  Fig.  10,  where  the  normalized  growth  rate  X-  is 

plotted  versus  t  for  R  u  /c=3.75,  R  „  / c= 2.4,  kc/w  =0.2,  y=1.118,  8  =0.4 

c  c  w  c  c  i 

and  v=0.002.  In  these  particular  parameters,  the  maximum  growth  rate 
X^=0.024  occurs  at  the  dielectric  constant  e=l.  The  range  of  e  corres¬ 
ponding  to  instability  occurs  repeatedly  as  the  value  of  dielectric  con¬ 
stant  e  is  increased  to  infinity.  For  example,  in  Fig.  10,  the  system 
is  unstable  for  the  dielectric  constant  e  satisfying  l<e<2.9  (first 
radial  harmonic  number  n=l)  and  10.2<e<13.2  (second  radial  harmonic  num¬ 
ber  n=2) .  In  this  value  of  conducting  radius  (R^cj^/c3! . 75) ,  without  a 
dielectric  material  the  system  can  be  unstable  only  for  first  radial 
harmonic  number  n=l.  Since  a  dielectric  material  reduces  significantly 
the  radial  wavelength  of  the  waveguide  mode,  it  is  possible  to  have  high 
radial  mode  perturbations  for  a  large  dielectric  constant. 

2 

Examining  carefully  on  the  cubic  equation  (46)  and  noting  (l+kc8z/wc)  > 
2 

(kc/w  )  in  a  typical  cyclotron  maser  instability,  we  conclude  from  Eq. 
c 

(46)  that  for  instability,  (b  +b  )  >0,  or  I (b  +b  )  |<<1  for  (b  +b  )  <0. 

o  o  o 

Moreover,  in  order  to  have  a  smooth  unstable  spectrum  with  a  broad  band¬ 
width,  it  is  required  for  the  ratio  a  to  satisfy 
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i  a 

for  a  considerable  range  of  k  space  when  the  beam  is  very  tenuous  (v<<0.01). 

For  present  purposes,  to  illustrate  the  influence  of  dielectric 

material  on  instability  bandwidth,  we  present  in  Fig.  11  plots  of  (a) 

(b  +b  )  ,  (b)  the  ratio  a  and  (c)  the  normalized  qrowth  rate  y.  versus 

+  w  1 

o 

kc/uj  obtained  from  Eqs.  (17),  (18),  (40)  and  (46)  for  R  w  /c=2.4,  y=1.118, 

c  w  c 

8  =0.4,  v=0.002  and  several  pairs  of  parameters  (R  w  /c ,  c ) .  Several 

X  c  c 

important  features  are  noteworthy  in  Fig.  11.  First,  the  matched  geo- 

metric  configuration  (R  w  /c,  t)=(3.75,  1)  is  most  unstable.  However, 

c  c 

the  range  of  k  space  corresponding  to  instability  is  narrow.  Second,  a 

slightly  unmatched  configuration  (R  to  /c,  e)=(3.75,  2)  is  unstable  for 

c  c 

a  broad  range  of  k  space,  with  a  considerably  reduced  growth  rate. 

Unfavorably,  the  growth  rate  for  (R  to  /c,  c)={3.75,  2)  is  a  sharply 

c  c 

changing  function  of  the  wavenumber  k.  Evidently  from  Fig.  11(b),  the 
ratio  a  for  (^  to^/c,  e)=(3.75,  2)  is  not  satisfying  Eq.  (47).  Third,  it 

is  noted  from  Fig.  11(b)  that  the  ratio  a  for  (R  to  /c,  <■. )  =  ( 3 ,  6)  satisfies 

c  c 

Eq.  (47),  thereby  resulting  in  a  smooth  unstable  spectrum  with  a  broad 

bandwidth  on  k  space.  Although  the  growth  rates  for  this  case  are 

slightly  reduced  from  those  values  of  (R  to  /c,  c)=(3.75,  1),  the  range  of 

c  c 

k  space  corresponding  to  instability  for  this  case  is  almost  twice  broader 

than  that  for  (R  w  /c,  t)=(3.75,  1).  In  this  regard,  the  unstable 
c  c 

spectrum  of  (R^oj^/c,  e)=(3,  6)  is  most  desirable  for  a  communicational 
application. 

We  conclude  this  section  by  emphasizing  that  stability  properties  of 
cyclotron  maser  instability  in  a  dielectric  loaded  waveguide  exhibit  a 
sensitive  dependence  on  the  dielectric  constant  e  and  thickness  of 
dielectric  material.  It  has  been  shown  that,  by  an  appropriate  choice  of 


b  +b 


w  —  (b  +b  ) 
c  dai  + 


■•0.1  , 


(47) 
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c  and  thickness  of  dielectric  material,  the  bandwidth  of  instability 
can  be  increased  more  than  twice  of  that  for  a  plain  conducting  wave¬ 
guide,  slightly  reducing  the  growth  rate  of  instability.  Moreover, 
dielectric  material  reduces  the  radial  wavelength  of  the  waveguide  mode, 
thereby  significantly  modifying  the  nature  of  parturbations. 
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VI.  CONCLUSIONS 


In  this  paper  we  have  investigated  the  influence  of  finite  wall 
impedance  effects  on  the  cyclotron  maser  instability  of  a  hollow  electron 
beam  in  a  waveguide  with  an  arbitrary  impedance  Z.  The  stability  analy¬ 
sis  has  been  carried  out  within  the  framework  of  the  linearized  Vlasov- 
Maxwell  equations,  assuming  that  the  beam  thickness  is  much  less  than  the 
equilibrium  radius  Rq  of  the  beam.  The  formal  stability  analysis  for 
azimuthally  symmetric  electromagnetic  perturbations  (3/36=0)  has  been 
carried  out  in  Sec.  II,  including  the  influence  of  finite  wall  impedance 
Z.  In  Sec.  Ill,  properties  of  the  vacuum  TE  mode  dispersion  relation 
have  been  investigated  in  the  absence  of  a  beam.  A  detailed  analytic  and 
numerical  investigation  of  cyclotron  maser  instability  has  been  carried 
out  in  Sec.  IV  for  a  wall  i  jedance  independent  of  the  eigenfrequency  and 
axial  wavenumber.  One  of  the  most  important  features  of  this  analysis  is 
that,  for  a  purely  resistive  wall,  the  growth  rate  is  substantially 
reduced  by  introducing  even  a  very  small  amount  of  resistivity  into  the 
wall.  Moreover,  the  range  of  kc/w^  corresponding  to  instability  increases 
rapidly  as  the  resistivity  of  wall  is  increased.  Cyclotron  maser  stability 
properties  in  a  dielectric  loaded  waveguide  have  been  investigated  in 
Sec.  V.  It  has  been  shown  that  stability  properties  of  cyclotron  maser 
instability  exhibit  a  sensitive  dependence  on  the  dielectric  constant  e 
and  thickness  of  dielectric  material.  By  an  appropriate  choice  of  e  and 
thickness  of  dielectric  material,  the  bandwidth  of  instability  can  be 
increased  more  than  twice  of  that  for  a  plain  conducting  waveguide. 
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Figure  lb  -  Contours  of  Constant  Phase  Angle  4>  and  Magnitude  |z|  of  Function  Z-Jj_  (xon> /xonJi_  <xon}  ln  the 
Plane  of  the  Root  x  =x  -x.  for  n=2 . 
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Figure  4  —  Plots  of  the  Normalized  Growth  Rate  Versus  kc/io  for  jz|=0.05,  Several  Different  Phase  Angles 
and  Parameters  Otherwise  Identical  to  Fig.  3.  C 


Parameters  Ident 


300  0  =  n  °  =  a/07  M 


/c  and  Parameters  Otherwise  Identical  to 


R^c  =  2.4.  R^c  =  au,y  =  1.118,  PL 
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Figure  lib  -  Plots  of  the  Ratio 


lRea jJc,  c)  =  (3.75.  1) 
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